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Abstract 

A theory of collective excitations in Bose-Einstein condensation in a trap is developed based on 
the quantum Hamilton-Jacobi equation of Bohm and the phase coherence along with the idea of 
off-diagonal long range order of Penrose and Onsager. First, we show that a free surface behaves 
like a normal fluid - a breakdown of superfluidity. Second, inside the free surface it is shown that 
the spectrum of phonons is of the form u = ck scaled with the external potential, where the speed 
of (first) sound, c = [A-Traph 2 ] 1 / 2 /M and k is the wave number. Third, in the limit a — > 0, the 
hard spheres in the Bose-Einstein condensation collapse to a close-packed classical lattice with the 
zero-point vibrational motion about fixed points. 

PACS numbers: 03.75.Fi, 03.65.-w, 03.65Sq, 05.30-d 
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The new Bose-Einstein condensation (BEC) in a trap is one of the remarkable phenomena 
exh.b.ted by a macroscopic system of Bose gas at low temperature Q. In his seminal paper 
on a theoretical model of a Bose gas |2j, Bogoliubov has laid out the basis for much of our 
understanding of superfluidity of a system of weakly interacting Bose particles. The concept 
of BEC is generalized to strongly interacting systems by Penrose and Onsager based on 
the notion of off-diagonal long range order (ODLRO) [3( which gives rise to superfluidity 
in He II. Many properties of an interacting Bose system can be understood in terms of 
independent quasi-particle excitations (phonons and rotons) of the system. These excitations 
are essential for a complete understanding of certain properties of the interacting Bose 
system, in particular, the longitudinal excitation spectrum in the phonon regime is crucial 
for a quantitative experimental confirmation of BEC in a trap. 

However, the Bogoliubov theory of collective excitations (first sound) differs so basically 
in its formulation from that of Feynman's atomic theory of two-fluid model ^| and Landau's 
two-fluid model |J that the problem of the collective excitations in BEC represents the major 
stumbling block in our attempt to confirm the new BEC in terms of the Bogoliubov theory. 
For instance, the new Bose condensates in a trap are extremely small in size (i.e., a liquid 
droplet of tens of microns in radius), and unlike He II they are highly compressible and have 
a nonuniform density profile. In addition, the new BEC is dilute (mean particle density p 
satisfies pa 3 -C 1) (7[. It is therefore obvious that the formal extension of Bogoliubov's theory 
to the new BEC is difficult on account of its inhomogeneous density profile. Furthermore, it 
is impossible in our study of the collective excitations to make use of the elegant quantum 
field theory technique developed in many-body theory which is valid only for a uniform 
system 

The new BEC is, however, the simplest of all condensed many-body entities from the 
microscopic point of view. The atoms in a trap may be treated as a weakly interacting 
Bose gas which can be accurately described by a self-consistent Hartree equation with the 
well-known scattering length in the hard sphere approximation for a pair- interaction jjj. 

On the experimental side, the first systematic experimental study on the collective ex- 
citation was recently carried out in a long 87 Rb Bose condensate in which the sound wave 
propagation can be treated approximately as a 1-D problem Yet the experimental data 
yield a qualitatively correct dispersion curve of the Bogoliubov spectrum u = ck with the 
estimated speed of first sound c e // = 2.0 ± 0.1mmsec~ l . Before the ideas of the new BEC 
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Q] in a trap become tenable, it will be necessary to show experimentally that the collective 
excitation spectrum in the new BEC indeed agrees well with the Bogoliubov spectrum with 
the proper geometrical corrections 

Theory, in the absence of credible experiments, could wander off down a blind alley |7[. 
And yet it is the theoretical insight that could unravel how Nature creates such a unique 
object with a peculiar property. Our analysis shows a remarkable new phenomenon - a 
breakdown of superfluidity at the free surface of BEC in a trap. This may be interpreted by 
broken gauge symmetry 10, 2, 3] and is precisely the reason why the speed of first sound 
was difficult to measure 9JI. 

n 

The object of this paper is, first, to show how Bohm's quantum theory can be applied 
to the problem of collective excitations in BEC in such a way as to incorporate Feynman's 
atomic theory of two-fluid model and the concept of phase coherence in ODLRO, and, 
secondly, to present a demonstration that symmetry breaking is taking place at the free 
surface with the first-order symmetry breaking perturbations in Lagrangian coordinates 
[tt! lis! ll^ . It is shown here for the first time that the dispersion relation for the surface 
wave is a manifestation of the broken symmetry at the free surface. The symmetry breaking 
also accompanies massless phonons (Nambu-Goldstone bosons) [13], the spectrum of which 
is shown to be the same in form as that of the Bogoliubov dispersion relation u = ck scaled 
by the external trapping force with a small geometrical correction, where the speed of (first) 
sound, c = lAiraph 2 ] 1 / 2 IM and k is the wave number [A. 

nn 

Since the perturbation techniques of quantum field theory jj, la] cannot be applied to 
a finite, spatially inhomogeneous system, we must find a new perturbation method that 
yields the Bogoliubov spectrum of collective excitations in the new BEC in a trap. Here 
we present straightforward, but mathematically precise, first-order perturbations by particle 
displacements in Lagrangian coordinates to an ensemble of particle trajectories defined by 
the condition of BEC that automatically include all the excited modes which tend to be 
otherwise omitted. In the present analysis, the problem of collective excitations reduces to a 
boundary-value problem for a BEC droplet in a trap bounded by a free surface. The present 
analysis also gives a good account of the possible formation of a close-packed lattice of a 
trapped Bose gas. 

The essential point of ODLRO is that BEC (and superfluidity) be described as a state 
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in which the density matrix can be factorized as 

p{x,x') = ip^(x)ip(x') +7O (1) 

where 7 — > as \x — x'\ — > 00. Here the single particle wave function ip{x) (mean field) 
represents the Bose condensed state [21] and j tp^ipdx = N, where N is the number of 
particles in BEC. With the the hard sphere approximation of the inter-particle repulsive force 
for Bose particles jsj], one can show that the mean field satisfies the nonlinear Schrodinger 
equation (Gross-Pitaevskii), 

Bib fr 2 

i% lH = ~2M V ^ + V{x)ext + (2) 

where g = A^h 2 a/M and a is the s-wave scattering length. Eq. is a self-consistent Hartree 
equation for the Bose condensed wave function. This short-range interaction brings about 
the Bose condensed state also ensures that the system possess the longitudinal collective 
excitations (phonons) in BEC. It should be noted that the nonlinear term \ip\ 2 is invariant 
under a U(l) group transformation. Hence we can apply Bohm's interpretation of quantum 
theory to Eq. (J2J) with the understanding that it is a part of the potential defined in his 
quantum theory. 

We write the equations for the ensemble average energy in the usual quantum theory 
00,0: 

n = / (-|^ v2 + V ( X U + f H 2 ) ^ ( 3a ) 
£ave = j (l^ 1 + V W««M 2 + f l^l 4 ) dx > j ^ dx = N - ^ 

At this point it is possible to show that the ground state density is given in terms of 
an external potential and the chemical potential upon minimizing the energy functional 
Eq. (j3b|) for a fixed number of particles in a trap with the condition p = (Penrose-Onsager 
criterion for BEC): 

p(x) = mx)\ 2 = -^ r [^-V(x) ext ], (4) 

where [l is a Lagrangian multiplier and is the chemical potential. It should be noted that the 
ground-state wave function Eq. (@J) has a nodal surface at the boundary where the symmetry 
of the Bose system breaks down. 
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To briefly recapitulate Bohm's interpretation of quantum theory [14], we write the wave 
function in the form ipjr.t) = f(r,t)exp[j:S(r,t)] = p 1 / 2 (r,t)exp{jrS(r,t)], where S(r,t) is 
an action (or a phase) |lj|. We then rewrite Eq. ((2)) to obtain, 

t + V-O^O ( 6 a) 

^ + (Z£)! + w a;) _il [ ZV_i(vrt! ] _ (5b) 

dt + 2M + { > 4M 1 p 2 J ~ ' 1 1 

where ^(cc) = V(x) ext + gp. The last term U eqmp of Eq. (}5b"|) is the effective quantum- 
mechanical potential (EQMP) which plays an important role in an inhomogeneous conden- 
sate and fluctuates near the free surface as will be shown shortly. 

From the point of view of macroscopic physics, the coordinates and momenta of an 
individual atoms are hidden variables, which in a macroscopic system manifest themselves 
only as statistical averages. In the following we show this interpretation by deriving the 
dynamical equations from Eqs. © with v(r ,t) = VS'(r ,t)/M for a single particle in 
condensation represented by the mean field in Eq. (j2J). 

In a dilute Bose gas in a trap, the force derived from the potential V(x) ex t balances 
the repulsive force due to gp(x) to create BEC for which p(x) = 0. Hence in the limit 
h — > or with the less stringent condition Vp = in a homogeneous medium for which 
U eqmp = 0, S(x,t) is a solution of the Hamilton- Jacobi equation in Bohm's theory 
Eq. (p)a| expresses a probability density of a statistical ensemble for a system. Eq. (j5| 
implies, however, that the particle moves under the action of the force which is not entirely 
derivable from the potential V(x), but which also obtains contributions from U(x) eqmp . It 
is the essence of Bohm's interpretation of quantum theory. In the classical description, the 
particle follows a classical orbit determined by the Hamilton- Jacobi equation. Eqs. © is a 
function of macroscopic dynamical variables and suggest us how we may choose the most 
stable particle trajectories (classical orbits) about which we may linearize the equations of 
motion. 

The first step in the development of a theory of collective excitations in BEC is to associate 
with a single particle wave function in ODLRO with precisely definable, symmetry breaking 
perturbations to the first order. Since external perturbations are always in the classical 
domain, the problem of collective excitations is reduced to finding atomic displacements 
that are consistent with the atomic theory of the two-fluid model of He II by Feynman in 
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which the exact partition function is given as an integral over particle trajectories, using 



a. 



We wish to stress that the collective 



his space-time approach to quantum mechanics 
excitation spectrum is also a test of the perturbation method for calculating the longitudinal 
excitations of (first) sound in BEC in a trap. 

Our perturbation method is based on the following physical picture: the collective coor- 
dinates are attached to the particles in motion as the phase coherent sound wave propagates 
from the center to the surface of BEC, and that the sound wave must also satisfy the 
prescribed initial-boundary conditions at the surface. The surface excitation involves the 
breaking up of phase coherence over the entire surface and is accompanied by a dissipation 
process of the sound wave. This gives rise a surface energy to conserve the total energy of 
an isolated system. 

To describe this, we introduce symmetry breaking perturbations to the particles in BEC in 
Lagrangian coordinates 17,QQ|, which can be defined as x = x + £(x , t). Here $,(x , t) 
is a function of the unperturbed position of a particle in the condensation, and remains 
attached to the particle in the condensate as it moves. This semiclassical perturbation 
method is in fact a well-defined accurate quantum mechanical approximation scheme which 
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yields the Bogoliubov spectrum for the collective excitations 2, 20J. 

At this stage we find it more convenient to work with the familiar dynamical equations 
instead of Eqs. (0). We also note that the use of the Hamilton- Jacobi equation in solving 
for the motion of a particle is only matter of convenience. 

If we write v(xo,t) = 'VS(xo,t)/M - (a solution of the Hamilton- Jacobi equation), we 
then obtain the following three dynamical equations from Eqs. and (jSJ): 

the equation of motion for a single particle, 

d 

M(—v + v-Vv) = -Vp, (6) 

the equation of continuity, 

|p + V • (pv) = 0, (7) 

where p is henceforth interpreted as the number density, 
and the equation of state, 

p(r, t) = pi oc [p(r, £)] + V ext , (8) 
where p\ oc [p(r, t)] is the chemical potential in the local density approximation j^ . 



Next we expand S, v, and p to the first-order, 

S(r,t) = S(r Q ,t)+£. V Q S(r ,t) (9a) 
v(x,t) = ^Z (9b) 
p(x, t) = p(x ) - V ■ \p(x )£] = p{x ) + 5p, (9c) 

where p(x ) = in BEC and Vo denotes the partial derivative with respect to x Q with 
V -> V - V £ • V . Eq. flUcJ) was derived by substituting Eq. (J9EJ to Eq. © with 
VS'(a;, t) = p(x, t) without taking into account of EQMP in Eq. (j5b|) . then integrating over 
time. 

Now it is easy to express from Eqs. (JUJ) the boundary conditions in terms of £ by the 
usual definition of incompressibility and irrotational motion of a fluid at the free surface in 
BEC j^, 0|; they are given mathematically by V • £ = and V x ( = 0. 

The phase coherence is essential to relate the Bose condensed wave function (mean field) 
in ODLRO to a many-body ground-state wave function, £ • VS , (x , t) = £\ ^ ■ 'ViSo^Xq^, t) 
in our analysis, where Sq^Xq^, t) is the phase of a single Bose particle in the system, and 
is a necessary and sufficient condition for phase coherence. Hence the collective excitations 
we study in this paper are phase coherent waves. 

The basis of this paper is the set of linear equations, Eqs. ©, ©, ©> an d Q that 
describe the collective excitations in BEC. It should be remarked that ip can support the 
collective excitations (phonons) only if it obeys the nonlinear Schrodinger equation for which 
we take 'VS(xo,t) = p(xo,t) as a solution of the Hamilton- Jacobi just as the nonlinear 
Maxwell equations describe radiation processes (photons). 

By virtue of Eq. (|9*c|) . the chemical potential may be expanded as 



d 

li(x,t) = /Uioc(cco) + V cxt - 7^Guioc)[V ■ (po£)]- (10) 

In BEC, p(xo) = 0, and the equation of motion gives 

Pioco(po(^o)) + V cxt = constant = pi oc o[po(0)] = p - (H) 

If we now set V ex t = u^r 2 /2, the density profile describes a spherically symmetric, nonuni- 
form Bose-Einstein condensation (SSNBEC). This is the problem from which our investiga- 
tion started. For the special example of collective excitations in BEC in trap, we limit our 
analysis to the problem of SSNBEC. 



In order to describe the collective excitations, we linearize Eqs. (J6J) along with Eq. (J9] 
Eq. (j^cj) . Eq. (fTUj) . and Eq. (fTT|) . It is straightforward algebra, though somewhat tedious, to 
arrive at the first-order equation of motion in The result is 

= M Va ~ ' V)r + (r ' V * )] ~ ^ [ar + r 2W] - (12) 
Here we have taken V ext (r) = MuJqT 2 /2, dfii oc /dp = 4irh 2 a/M, a — V • £, and have also 
dropped the subscript in r . 

Eq. (fT2"j) is a typical second-order (inhomogeneous) partial differential equation subject to 
the initial-boundary conditions and has two particular solutions, one corresponding to the 
surface waves and the other to longitudinal sound waves. The general solution to Eq. (J12j) 
is a combination of the two solutions. 

We begin with the surface phenomena, and shall study later the more complex problem 
of the longitudinal sound waves in SSNBEC. We limit the surface waves by imposing the 
boundary conditions V • £ = and V x £ = 0, which simplifies the algebra considerably. 
This implies that we may solve for the potential flow as V 2 x = and £ s = — V%. Here the 
subscript s stands for the surface waves. The solution for x is given by 

X(r,t) = Y^[xi(ty + X i(t)r-^]Y l>m (0,<t>), (13) 

£,m 



where we may set x~ (t) = for a quantum liquid 



of three orthogonal vector spherical harmonics 



droplet. We may now expand £ s in terms 



24] 



& = EtdrMK+^rM)^ +& m Ma 3 ]. (u) 

Here we have defined the three vector spherical harmonics as 

ax = e r Y Am (0,</>), a 2 = rVY £ym (8, 0), a 3 = r x V^, m (0». 

It is now only a matter of elementary algebra, by substituting Eq. (|T4|) into Eq. (|T2|) to 
obtain a dispersion relation, 

^L/ = ^o- (15) 

Here we have taken x+{t) = elwt - It is a ^ once apparent that the dispersion relation for 
the surface waves is independent of the internal dynamics of the trapped Bose gas (e.g., the 
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pair interaction potential). Another way of looking at the property of the free surface is 
to note that the surface waves are driven by the action of external perturbations over the 
surface of a spherical droplet in equilibrium, similar to that of gravity waves on the surface 



of a fluid in equilibrium in a gravitational field [25|, |26[. More importantly, the dispersion 
relation is independent of H which implies that the free surface behaves like a classical fluid. 

It is, therefore, evident that Eq. (fT3j) is the first definite proof that the free surface of a 
superfluid behaves like a normal fluid; it is a peculiarly universal property of the free surface 
of a superfluid under an external field. Furthermore, Eq. ()15|) is also an excellent example of 

Hf . The breakdown of superfluidity has been already observed 
23, Q|- The breakdown of superfluidity at the free surface 



broken symmetry in Nature 
at the vortex core in He II 

in rotating He II (^(J also resolves a long standing puzzle in Landau's two fluid model in 
connection with the question of why a superfluid component rotates with the normal fluid 
at the free surface p, |29j, |30( • 

Since in practice the free surface is never confined to a strictly mathematical surface, the 
thin surface layer which also contains phonons and particles driven by EQMP, 



U = .«L|!V.iM ]= .m ( 16 ) 

Thus the particles experience a force from U eqmp which fluctuates with the momentum 
p = ViS and energy of a particle near the surface with the degree of divergence Muj 2 fj, /D 2 
as D = \p — (l/2)Mu 2 r 2 ] — > at the free surface as emphasized by Bohm 14]. U eqmp , 



therefore, breaks up phase coherence in the collective excitations in the surface layer. 

In our model of an imperfect Bose gas, between each pair of particles there is a hard sphere 
repulsion of range a and no long-range interaction. This pair-interaction ensures that the 
system possesses the longitudinal collective excitations (phonons) in BEC. Thus the problem 
of collective excitations in SSNBEC should be considered as a spherical longitudinal, phase- 
coherent sound wave propagation. In the phonon regime, the dispersion relation for phonons 
must be of the form u = ck scaled with the external trapping force, where the speed of first 
sound, c = [Aiiaph 2 ] 1 / 2 /M, and k is the wave number. Because of our approximation for 
the pair-interactions of Bose particles in the system, the presence of the s-wave scattering 
length a in the speed of (first) sound wave c is essential. Also this approximation should 
be valid so long as the condition pa 3 <C 1 is met, a condition which restricts the discussion 
to the region of temperature near absolute zero, which explains why there cannot be roton 



excitations in BEC in a trap. 

Now we are ready to define the compressibility of the fluid as o = V ■ £ from Eq. (|5cj) . 
and then its spectrum can be obtained from Eq. (J 12)) by taking the divergence on both sides. 
After straightforward algebra with the help of the vector identity V • [(r ■ V)£] = r ■ Vct + a 
and the condition of superfluidity V x £ = 0, we obtain 

§p*(r, t) = l -u 2 {a 2 - r 2 ) VV - 3uj 2 r-^a - 5u 2 a, (17) 

where a 2 = 8irh 2 ap (0)/(Muj ) 2 . 

The solution of the equation is not entirely trivial. Writing a(r,t)= S(t)W(r)Yi jTn (8, 0), 
we obtain the variable separated equations: 

^S(t) + X n S(t) = (18) 

and 

(« 2 - **)[~(rW n (r)) - ^±^-W n (r)] 

d 

-6r—W n {r) 

+ (2X n /uj 2 -10)W n (r) = 0, (19) 

where A n is a constant of separation. The eigenvalues are determined by Eq. ()19j) . and by the 
boundary conditions on a(r, t), i.e., cr(r,t) = at the free surface and the origin. Here we 
tacitly assume an infinitely small point source at the origin to drive an outgoing spherical 
sound wave [25[. The speed of first sound c = [Anap^Ti 2 ] 1 / 2 / M is the maximum at the 
origin and becomes zero at the free surface. To see this is indeed the case we recall that the 
surface layer is no longer a superfluid but is a normal fluid. Hence the phonons can interact 
with the normal fluid and, upon the dissipation process, give rise to a surface energy. 

The solution of Eq. (|18J) then gives the dispersion relation for the collective excitations 
with the eigenvalues that are a function of the s-wave scattering length, the radial trap 
frequency, the peak density at the center of the trap, i.e., the speed of first sound c ctr = 
[47rap(0)fr 2 ] 1//2 /M, the trapping frequency, and the wave number with kg ~ Ijr. To show this 
explicitly, we transform Eq. ()19j) to the Sturm-Liouville problem and obtain the eigenvalues 
in terms of a complete set of functions with the orthogonal properties: 

An = [ b r 2 (a 2 - r 2 f[{^W n ) 2 + l -^^W 2 n {r)}dr (20) 
Jo dr r 
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and with 

r-6 



r 2 (a 2 - r 2 ) 2 W m (r)W n (r)dr = 6 m , n . (21) 

Here A n = (2X n /uj 2 — 10) and b = a = [8nah 2 p^)] 1 ^ 2 / (Muq) is the radius of the condensate. 
This procedure, though somewhat a digression, is the only way we can explicitly show the 
functional relation between the eigenvalues and the speed of first sound. 

Equations (|2()jl -(|21 jl exhibit the main feature of our method: that the eigenvalues are, 
indeed, a function of the speed of (first) sound, b = [\/2/loq\ c ctr , where the sound speed c ctr = 
[47rap(0)/i 2 ] 1//2 /M, and the trapping frequency too, for all values of the angular momentum 
i. This method of presentation also shows why the excitation spectrum has not shown the 
dependence, or lack thereof, the speed of (first) sound in previous studies 7 1 ]. The functional 
relation of the speed of (first) sound and the wave number in the excitation spectrum is also 
clear in this integral representation. More importantl y it is apparent that the perturbation 
method by Lagrangian displacement vectors 

HQ 

19] is in fact a well-defined quantum 

mechanical approximation scheme. 

Another important point is that the eigenvalue is a function of the speed of first sound not 
at the free surface but at the center of SSNBEC. This has a simple physical interpretation: 
an outgoing spherical sound wave can travel with little energy loss toward the free surface 
where the sound wave (phonons) interacts only with the normal fluid of the surface layer 
and completely dissipates giving rise to a surface energy. It is precisely the nature of a 
superfluid that it cannot interact with phonons, but the phonons will interact with a normal 
fluid as Mott emphasized [31]. Therefore, the phenomenon of the fluctuation of particles 
due to U eqmp and the dissipation of sound waves at the surface layer can be understood in 
terms of Kubo's fluctuation-dissipation theory 0] and satisfies the conservation of energy 
for an isolated system. We thus see that the particles in the surface layer obey the Boltzman 
statistics just like rotons because of the elevated particle energy. Hence the surface layer is 
a normal fluid in every respect jsj]. 

The analytical expression Eq. (j20J), in principle, gives a complete solution to the eigenvalue 
problem. A moment's thought shows, however, that the expression for W n is too complicated 
to evaluate the eigenvalues by this approach. Here we follow a different, but equivalent, 
approach to the eigenvalue problem; that is, we shall solve Eq. (fTTIj) as an eigenvalue equation. 
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With the substitution W n (r) = r^ i+1 ' 2 ' l / 2 Z n (r) together with x = r 2 /a 2 , we obtain 

d 2 d 
x(l - x) — Z n + [c - (a + b + l)x\—Z n - abZ n = (22) 

This is just the Gauss differential equation [33] with c = ±(£+1/2) + 1, a+b = ±(£+l/2) + 3, 
and ab = (1/4) A n - 6 [±(£ + 1/2) - 1/2]. 

Eq. (|2*2*jl has regular singular points at x = 0, x — 1 and x = oo. Its solution is the 
hypergeometric function, which is analytic in the complex plane with a cut from 1 to oo 
along the real axis j^ . 

Since we are interested in the low-lying excited states (phonons), it is only necessary to 
find the smallest eigenvalues in the differential equation. This is consistent with Feynman's 
picture of phonons - a sound (longitudinal) wave without nodes |4]. A simple, but accurate, 



numerical method 34| has been employed to evaluate the eigenvalues in the domain [0, 1] in 
which the solutions are analytic. 

Now returning to Eq. |T8j) and taking S(t) = e luJt , we obtain the dispersion relation as 

u ph = ±[A s /2] 1 / 2 w , (23) 

where X s is the smallest eigenvalues from Eq. f!22D . 

The principal result of this paper is Eq. (|23|) . It is given in Fig.l, where the ratio u p h/uJo 
with respect to the angular momentum i is plotted. A close examination of the disper- 
sion curve shows that in the phonon regime the energy spectrum of longitudinal collective 
excitations is nearly linear with respect to i (i.e., the wave number kg ~ ijr). 

In an inhomogeneous medium, the Bogoliubov dispersion relation u p h = ck with 
c = [47rap(r)/t 2 ] 1 / 2 /M must be studied approximately with fce ~ f/r at a given radius 
in SSNBEC. In particular u p h/ujQ = 1.5972 for £ = is a unique value in a finite space 
problem, which corresponds to a uniform radial perturbation. Finally, since the eigenvalues 
are positive for all values of £, the Landau critical velocity v c = e p h/p is also finite. Hence 
the imperfect Bose gas is a superfluid, whereas an ideal Bose gas, for which c is independent 
of the scattering length, is not, because c and v c vanish identically. Moreover, the dispersion 
relation, Fig.l, also shows how the presence of the repulsive pair-interaction determines the 
excitation spectrum of phonons, which is nearly linear with respect to kg 

As we have shown above, the symmetry of the ground state wave function breaks down 
at the free surface, i.e., a nodal surface. It is therefore natural to identify the underlying 
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basic mechanism for the symmetry breaking as a spontaneously broken gauge symmetry at 
the free surface which accompanies phonons as Nambu-Goldstone bosons |l2|, [l3| . 

Next we discuss a possible formation of a close-packed lattice as a — > in BEC. Our 
ground state wave function describes a perfect Bose gas modified as little as possible by the 
presence of the repulsive pair-interactions (a hard sphere approximation). Hence the hard 
spheres in BEC could collapse to a close-packed lattice in the limit, {i.e., a — > 0) in which 
each particle vibrates about a fixed center. Since the number of particles N in the system 
must be conserved in the limit, the density must satisfy (4/3)7TG^p c z = N; p d is the mean 
density of the collapsed spheres 0,13- ft i s unlikely that the closed-packed lattice will be 
realized in a trap, because the radius of the condensate b — > d c i <C 1 as a — > [Eq. (|2(J|) ]. 
which implies a far stronger external trapping force to overcome the inter-particle repulsive 
force. However our analysis may provide a qualitative mechanism by which the close-packed 
lattice can be formed from BEC by symmetry breaking. 

The separation of the zero-point vibration is not always possible, but in the limit a — > 0, 
b = a = [87rah 2 p(O)] 1 / 2 / \Mujq) — > d c i, X n ~ 0, we have the following zero-point vibrational 
motion about fixed points. Since A n = 5.0uq, we have a dispersion relation for the closed- 
packed lattice from Eq. (JTSj) . 

wiatt = (5) 1/2 ^ , (24) 

which is again independent of h. This implies that if there is actually a transition from 
BEC to a crystalline ground state, it is a classical lattice. It is gratifying to note that the 
condition under which the close-packed lattice could be realized in BEC is essentially the 
same as that of a crystalline lattice observed under an external high pressure applied to He II 
The physical mechanism by which the close-packed lattice can be realized is again 
the spontaneously broken gauge symmetry which accompanies phonons as Nambu-Goldstone 



bosons 



12, 



In summary, we have shown the perturbation method with the Lagrangian displacement 
vectors on particle trajectories in the many-system yields the results that are consistent with 
those of quantum field theory technique of Bogoliubov. The compressible perturbations give 
the dispersion relation for a spherical (longitudinal) sound wave which is of the form of the 
Bogoliubov spectrum for phonons. The incompressible perturbations drive the free surface to 
behave like a normal fluid. The dispersion relations Eqs. (j!5|) and ()23|) are interpreted as the 
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result of the spontaneously broken gauge symmetry at the free surface which accompanies 
phonons (Nambu-Goldstone bosons). The present approach gives a unified picture of the 
collective excitations in BEC in a trap and also suggests a possible mechanism by which 
BEC reduces to a close-packed lattice under external pressure similar to that of the helium 
crystals at constant pressure. 

Note added in proof. After this paper was written, I learned about 't Hooft's conjecture 
that is essentially equivalent to our picture of fluctuation-dissipation of the longitudinal 
sound waves at the free surface, which is driven by the effective quantum potential and is 
due to the broken symmetry [Massimo Blasone, Petr Jizba, Giuseppe Vitiello, Phys. Lett. 
A 287 (2001) 205]. 
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Angular Momentum 



FIG. 1: The ratio ujph/ooo = [2A,,] 1 / 2 is plotted against the angular momentum i. It shows how the 
energy spectrum of phonons varies with the angular momentum. 
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